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With anomalies found in cosmic ray observations and unsolved questions of the Standard Model of 
particle physics like the discrepancy in the muon's anomalous magnetic moment, the idea of an (7(1) 
extension of the Standard Model arose. This extension allows for the interaction of Dark Matter by 
exchange of a photon-like massive force carrier A' not included in the Standard Model. We discuss 
the possibility to constrain the dark photon parameter space by using data taken from rare kaon 
decays. Therefore we analyze the decay K — > fivA' as a signal process and calculate the expected 
Standard Model background. Using these results we calculate new limits on the A' parameter space, 
providing motivations for new rare kaon decay experiments to extend the existing A' bounds. 
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The possibility to explain anomalies in astrophysics 
and particle physics [JIH by extending the Standard 
Model of Particle Physics (SM) by an additional U(l) D 
gauge group manifesting itself in a massive gauge boson 
A' ("Dark Photon") in the MeV to GeV mass range 0- 
6] motivated a strong activity in theoretical as well as 
in experimental physics @-[l6|. In a recent work [l5| . 
where new physics scenarios were explored as a possi- 
ble explanation of the proton charge radius problem, the 
constraints from rare kaon decays were examined. In the 
present work, we will use these rare kaon decays to find 
constraints on the A' parameters. The A' couples to the 
electromagnetic current via kinetic mixing giving rise to 
a QED-like vertex term 

where e is the kinetic mixing parameter describing the 
coupling strength to the electromagnetic current by 
e 2 = a' D /aQED with qqed = e 2 /4ir [9|. 
In the following we will study the process K + — > y&v^A' 
as a possible signal from the dark sector (see Feynman 
diagram in Fig. [I} within the mentioned framework 
of kinetic mixing (model I) as well as in a model 
where the A' couples only to the muon assuming an 
explicit breaking of gauge invariance applied in Ref. 
(model II). In the pioneering experiment [17| of the 
decay K + — > /i + + neutrals, only the charged muon is 
detected, excluding further charged particles or photons 
in the final state. Therefore this process cannot be 
distinguished from the K + — > y^v^vivi decay (Feynman 
diagrams in Fig. [2j , making it crucial to obtain a precise 
understanding of the background. This analysis will 
allow us to compute bounds on the A' coupling strength 
e as a function of its mass parameter ttia 1 ■ 

In calculating the decays of Figs. [TJ [5] we will use, that 
the largest external momentum scale of the considered 
processes, the kaon mass ttlk , is far below the weak gauge 
boson masses mw and mz allowing to approximate the 
weak interaction processes by point processes, respec- 
tively. Furthermore we parametrize the hadronic weak 
K + decay current by (fj|?vy a (l - j 5 )u u \K + (k)) = fxk a , 
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Figure 1. Feynman diagram for the process K — > ^v^A' 



where fx is the kaon decay constant. 
Since the K + is assumed to decay in rest the analysis 
will be performed in the K + rest frame where the kaon 
four- vector is given as k — (m^, 0). Furthermore we 
use the momentum four- vectors I — (E^, 0, 0, \l |) for the 
muon, q = (E V ,E V sin 9„, 0, E v cos9„) for the neutrino, 
and q' = k — I — q for the A', where we choose the z- 
axis along the muon 3-momentum. All other 3-momenta 
and all angles associated with other particle are not ob- 
servable and therefore these decays have isotropic muon 
angular distributions. 

The A' production amplitude derived from diagram[T](a) 
is given by 
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where Gf is the Fermi constant and 9 C the Cabbibo 
mixing angle. For model II the signal process is com- 
pletely described by this amplitude. Due to the need for 
gauge invariance within the kinetic mixing framework, in 
model I the amplitude is given by the coherent sum over 
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Figure 2. SM background from K + — > v^vivi. 



all diagrams of Fig. Q] 
Mm = 
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whereas W — k — q' and e 0123 = 1. The term propor- 
tional to fx is known as inner Bremsstrahlung contribu- 
tion (IB) and does not contain any structure effects. The 
contribution proportional to H pu contains the structure 
dependent terms, which are parametrized by the form 
factors Vi, A\, and A?.. 

In the K + rest frame the differential decay width for 
K + — > fi+i/pA' using the conventions of Ref. [HI then 
reads as 



dT(K+ -> n+v^A 1 ) 
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where we have fixed the angle V by evaluating the energy 
conserving (^-function as: 
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To obtain the decay width, Eq. ([3]) has to be integrated 
over E^ and E v within the limits as given in section IV 
of Ref. flBI . The background contribution in terms of its 
5-fold differential decay rate is given by 
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in the kaon rest frame, where 
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is fixed by the energy conserving <5-function, E 2 — p 2 and 
M is the coherent sum over the amplitudes derived from 
the Feynman diagrams in Fig. [2] 

In order to obtain a dimensionless quantity, it is helpful 
to consider the ratio of these decay rates relative to the 
one for the K + — > fJ. + v^ expressed by 
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Similar to the distributions ((3]) and ((JJ) this decay rate is 
isotropic in the muon angles. Thus all corrections due to 
an angular detector acceptance being not 47r will cancel 
each other in the ratio with respect to the decay rate 
T(K+ -> n+v„). 

For our analysis, we use the existing data published by 
Pang et al. in Ref. |17| . who found an upper limit of 
T(K+ fj+X)/T(K+ ->■ ijl+v^) < 2 • 1CT 6 on the ratio 
of the decay rates, where X is a possible, not detectable 
neutral state not containing a photon. In order to obtain 
the experimental limits from these data the differential 
decay rate H + v^A') has to be folded with the 

detector efficiency D(E fl ) given in Ref. [17j |. i.e. 
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Since the kinetic mixing factor e is a global fac- 
tor of the amplitudes ([1]) and ([2]), one can rewrite 
R (wia') — e 2 R (niA') and thus finds an upper bound 
for allowed values of e 2 as: 



e 2 < 



2 ■ 10~ 6 
R{m A ') ' 



(6) 



In Fig. [3J (upper and middle panel) the differential 
decay rate for the signal process relative to the decay 
K + — > is shown calculated within model I and II 

for the full phase space (left panels) and with applied 
corrections due to the given detector acceptance (right 
panels) , according to the experimental set-up of Ref. [17| . 
One notices that within the kinetic mixing model (upper 
panels of Fig. |3J the internal bremsstrahlung contribu- 
tion completely dominates the result for the considered 
A' mass parameters: comparison between IB curves and 
curves including the form factor dependence which was 
evaluated according to Refs. (lil. Il9|. Since in model II 
the gauge invariance is not required, the decay rate is 
enhanced by a factor of l/m A , compared to model I. The 
expected SM background from the decay K + — > n + v^vv 
with the applied experimental cuts of Ref. (l7| is shown 
in the lower left panel. 
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Figure 3. Upper and middle panels: Ratio of -^-(K + — > [i^v^A') and Y{K + — > p^v^) for various A! masses for perfect 

, , 

detector efficiency (left panels) and for finite detector efficiency of Ref. 17] (right panels) at e 2 = 1. Upper panels: kinetic 
mixing model (model I); middle panels: model II, where the A' only couples to the fj, + . Lower left panel: Standard Model 
background for different neutrino families using the detector efficiency function. Lower right panel: ratio of total decay rates 
T(K + fi+v^A') at e 2 = 1 relative to T(K + ->• v+v^vv) in model II. 



As one can see from the lower right panel of Fig. [3J the 
total A' decay rate (model II) T(K + — > fiv^A') calcu- 
lated with e 2 — 1 is about a factor of 10 9 larger than 
the decay rate to SM particles T(K + — » [i^v^vivi). This 
corresponds to an A! signal, which will dominate over the 
expected SM signal for mixing factors down to e 2 ~ 10~ 9 . 



The calculated limits on the A' parameter space 
are shown in Fig. 2] In this figure the colored regions 
correspond to already excluded configurations of mass 
and coupling strength 0, [l2|, EH ■ In this plot we have 
included the old Q as well as the new 2^, 21 1 exclusion 
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Figure 4. Exclusion limits on the A' parameter space. Dashed-dotted curve: bound calculated in the kinetic mixing model 
(model I) for an accuracy of the ratio T(K + -> A')/T(K+ of 10~ 9 . Dashed curve: result for the 1973 data |l7l | 

within model II, where the A' only couples to the fi + . Dotted curve: bound calculated in model II for an assumed improvement 
of the experimental accuracy by two orders of magnitude, i.e. 2 • 10 -8 . 



limits from (g — 2) of the electron compared to the fine 
structure constant a. The red colored region represents 
the so-called (g — 2) M welcome band, where the A' 
contribution to (g — 2) M could be invoked to explain the 
existing discrepancy |7fl. 

A possible bound for the kinetic mixing model 
is represented by the dash-dotted curve for 
an assumed experimental accuracy of the ratio 



A')/T{K^ 



H+u^) of 10- 



Based 



T(K+ 

on the old (g — 2) e exclusion limit, the 1973 data |l7| 
allow to slightly improve the bound at low masses and 
large e (dashed curve) within model II. Due to the 
refinement of the theoretical determination of (g — 2) e 
the bound from rare kaon data is already covered by 
the new (g — 2) e limit. Furthermore in Fig. [5] we give 
an estimate in which way the exclusion limits change 
due to an improvement in the experimental accuracy 
of the ratio T{K+ [t+v ll A')/T(K+ -> fj, + v M ) by 
two orders of magnitude (dotted curve). Obviously an 
improvement of the experimental quantities on the right 
side of Eq. ^ will allow the exclusion of a large region 
of the parameter space up to masses of about 80 MeV 
since the bound on e 2 is depending linearly on the RHS 
of Eq. © . Larger angular and momentum acceptancies 
and a larger rate of stopped K + compared to 17J for 
example will allow to improve this quantity significantly. 
Such an improved extraction might be achieved by new 
facilities, such as the NA62 experiment at CERN or rare 
kaon decay experiments at JPARC. 



We have used the rare kaon decay K + — > [i^v^A 1 to 
find a bound on the A' parameters extending the ex- 
cluded region at low masses in a model where the A' cou- 
ples only to the muon. We have shown, that the method 
used in this work may be suited to extend the existing 
limits within two models for the A' coupling. For that 
purpose more precise data are necessary. Improving the 
accuracy compared to that of the forty year old pioneer- 
ing work of Ref. [l7T | by two ore more orders of magnitude 
would allow the exclusion of a significantly larger part of 
the up to now allowed parameter space, which is also 
containing a considerable part of the (<? — 2) M welcome 
band. 
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